Charge pumping through the interface between topological insulator with massive 
surface Dirac fermions and a precessing ferromagnet 
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We study charge pumping through the interface between the ferromagnetic (F) overlayer and 
a three-dimensional topological insulator (TI) using multiphoton nonequilibrium Green function 
formalism. The microwave- driven magnetization precession of the F overlayer leads to the time- 
dependent exchange field which can induce surface energy gap via the magnetic proximity effect, 
thereby forcing Dirac fermions on the TI surface to acquire a time-dependent mass. This leads to 
the huge enhancement of the dc pumping voltage V^pump, which is proportional to the microwave 
frequency uj^ when compared with the case without the surface gap opening where Vpump is governed 
only by the spin-orbit coupling on the surface of TI. The predicted effect survives even when the 
Fermi energy is outside of the surface gap or with finite density of charge carriers in the bulk of the 
TI slab. 
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The recent experimental discovery [1-3 of three- 
dimensional (3D) topological insulators (TIs), which 
posses a usual bandgap in the bulk while supporting ex- 
otic metallic states on their surfaces, has opened new av- 
enues to explore low energy quantum phenomena in two- 
dimensional condensed matter systems, as well as bench- 
top setting to create particles (such as axions, monopoles 
and Majorana fermions [2]) conjectured in high energy 
physics. The strong spin-orbit coupling (SOC), as the key 
ingredient in this novel materials, causes spin- momentum 
locking of low-energy quasiparticles described by a sin- 
gle Dirac cone on the TI surface which is not seen in 
related systems, such as two-dimensional electron gases 
(2DEGs) [4 or graphene [5^ with the Rashba SOC. 

Since helical surface states of 3D TIs are gapless, they 
respond strongly to any time-reversal symmetry breaking 
perturbations such as covering the surface by magnetic 
impurities [6] or a ferromagnetic (F) overlayer [7 . In 
such situation, the Dirac cone is perturbed by the prox- 
imity induced exchange coupling which leads to a surface 
energy gap Agurf at the Dirac point observed in recent ex- 
periments [6 . The corresponding massive Dirac fermions 
exhibit half-quantized Hall conductivity, gh = ±e^/2/i, 
even in the limit of vanishingly small Agurf ^ 0, as long 
as the Fermi energy lies within the surface gap Agurf in 
clean samples [2 (or within the bulk gap in TI sam- 
ples where magnetic disorder localizes surface quantum 
states [8^). This effect, known in high energy physics as 
parity anomaly, is determined by the bulk topology and, 
therefore, serves to distinguish TI surface states from all 
other 2D systems or topologically trivial surfaces in con- 
densed matter physics [2 . 

However, unlike the well-known integer quantum Hall 
effect, where chiral edge states carry current and connect 
to external leads in mult iter minal Hall bars, cfh on the 
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FIG. 1: (Color online) Schematic view of the topological 
insulator-based semi-MTJ which pumps charge current oc uo 
along the x-axis due to magnetization dynamics driven by 
the absorption of microwaves of frequency uj under the FMR 
conditions. The semi-MTJ contains a single F layer of finite 
thickness whose magnetization m(t) is precessing around the 
z-axis. The precessing F layer is attached to the left semi- 
infinite N lead and to 3D TI slab on the right, so that as 
soon as the precession cone angle becomes finite energy gap 
Asurf(t) can appear on the surface of TI due to proximity 
effect governed by the component of the magnetization per- 
pendicular to the F|TI interface. Each layer is composed of 
monolayers (MLs) modeled on an infinite square tight-binding 
lattice. No dc bias voltage is applied between the N leads. 



surface of TI cannot he measured by conventional multi- 
terminal dc transport setups since either edge states do 
not exist (if the whole surface of TI is covered by mag- 
netic impurities) or one finds the usual cjh = /h when 
chiral gapless edge state propagates along the domain 
wall if only portion of the TI surface is covered by mag- 
netic impurities [2 . Thus, experimental efforts to detect 
manifestations of gapped helical surface states of 3D TI 
have been directed toward magnetoelectric effects probed 
using electromagnetic (EM) radiation in order to measure 
magneto-optical Kerr and Faraday rotation angles 
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for EM wave propagating through a slab of TI that is 
weakly exchange-coupled to F layer [9HTT]. Demonstrat- 
ing that for TI slabs thinner than the radiation wave- 
length is quantized in units of the fine structure con- 
stant while Ok = —90° would confirm the topological 
property of the bulk [11 . However, these experiments 
require negligible bulk conductivity — otherwise the uni- 
versality of and 6k values is broken [11^. This is dif- 
ficult to achieve in presently available TI samples which 
are always unintentionally doped (e.g., Bi2Si3 is n-doped 
by the Se vacancies), thereby obscuring [T2^ topological 
properties anticipated for perfectly insulating bulk [2 . 

In this Letter, we predict a novel charge pumping phe- 
nomenon in vertical F|TI multilayers in the absence of 
any dc bias voltage which is very sensitive to the pres- 
ence or absence of the surface gap Ag^rf and, therefore, 
can be employed to detect the presence of massive Dirac 
fermions by measuring the dc pumping voltage. Our de- 
vice setup is illustrated in Fig. [l] where finite thickness 
F layer in contact with finite thickness TI slab are sand- 
wiched between two normal metal (N) leads to form a 
semi- magnetic tunnel junction (MTJ). The dynamics of 
the magnetization of F layer as coherent macrospin pre- 
cession is induced by the absorption of microwaves under 
the conditions of ferromagnetic resonance (FMR). Such 
FMR situation in conventional ferromagnetic multilayers 
leads to pumping of pure spin current by the F|N in- 
terfaces, which has been amply explored over the past 
decade as a robust and ubiquitous effect at room tem- 
perature p!3l . 

On the other hand, using recently developed [14] 
nonequilibrium Green function (NEGF) theory of spin 
pumping in the presence of strong inter facial SOC, we 
demonstrate in Fig. [2ja) and (b) that proposed device in 
Fig. [l] would pump time-dependent charge current /(t), 
proportional to the microwave frequency cj, in the di- 
rection perpendicular to the interface F|TI. As soon as 
the cone angle becomes non-zero due to microwave ab- 
sorption, the time-dependent exchange field acquires a 
component {A sin sin ujt)ex which is perpendicular to 
the surface of the TI and can induce the corresponding 
surface gap Asurf(^). Note that the value of A^^^ is 
not directly related to A|sin^| since magnetic proxim- 
ity effect can be influenced by the properties of the F|TI 
interface (discussed below in relation to Fig. [3|. 

The detection of parity anomaly on the (non- 
localized [8^ ) surface of TI requires to measure a physical 
quantity which experiences a sudden jump in its value at 
arbitrary small Agurf and for the Fermi energy within the 
surface gap \Ef — £^dp| < Agurf (^dp is the energy of 
the Dirac point). Although we do not see such disconti- 
nuity in the setup of Fig. [l] we do find in Fig. |2|a) that 
dc pumping voltage in an open circuit Ipump = I/G{6) 
is about three orders of magnitude larger at small cone 
angles if A^^^ ^ when compared to the junction 
with Agurf = 0. Here / is the time- averaged value of I{t) 
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FIG. 2: (Color online) The dc pumping voltage in F|TI|N 
junction in Fig. [l] without Agurf = and with A™^? = 
0. 2 1 sin ^ I eV energy gap on the surface of TI when Fermi 
energy is: (a) inside the surface gap Ep — ^dp = 0.01 eV (for 
> 6°); or (b) outside of the surface gap Ef — E^p = 0.35 
eV for all 0. Panel (c) shows tunneling anisotropic magne- 
toresistance (TAMR) corresponding to situation in panel (a), 
while panel (d) plots dc pumping voltage as a function of the 
Fermi energy for precession cone angle 6 — 10°. 



in the N leads, and G{0) is the conductance of F|TI|N 
semi-MTJ computed by tilting the static magnetization 
of the F layer by an angle away from the z-axis and by 
applying the linear-response bias voltage between the N 
leads. The precession cone angle can be controlled exper- 
imentally via the input microwave power in the typical 
range < 30° [15], and F|TI heterostructures have al- 
ready been fabricated , so that confirmation of the re- 
sult in Fig.|2]is within the reach of present experimental 
technologies. Figures[2|b) and [2|d) shows that the signal 
l/pump(A^^f ^ 0) remains larger than Vpump (Agurf = 0) 
even when \Ef — £^dp| > A^^^, but pumping voltage in 
this case for the experimentally relevant range < 30° is 
much smaller than in the opposite case [except very close 
to £^DP where it vanishes in Fig. [2|d)]. 

The non-zero pumping voltage at the adiabatic level 
cx uj injunctions with interfacial SOC [14 (in the absence 
of interfacial SOC, simple F|I|N junctions pump [16] 
charge current oc cj^ which is, therefore, nonadiabatic) is 
closely related to the tunneling anisotropic magnetoresis- 
tance (TAMR). The out-of-plane TAMR for F|TI|N semi- 
MTJ is defined [14 as TAMR (0) = [i?(90° - 0) - R{e = 
90°)]/ R{0 = 90°) using conventional resistances mea- 
sured by tilting the static magnetization of the F layer 
in Fig. [T] away from the transport direction = 90°. The 
TAMR curves plotted in Fig. [2|c) show that change in 
the conductance G{0) = 1/R{0) is too small to account 
for the enormous difference between ypump(A^^£ ^ 0) 
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FIG. 3: (Color online) The local density of states 
N{E^ ky,kz = 0) on the surface of TI slab of thickness dp = 10 
MLs, as well as on the second and third ML marked in Fig.[l] 
which is either (a) weakly (7c = 0.15 eV) or (b) strongly 
(7c = 0.4 eV) coupled to the neighboring F over layer of thick- 
ness c^F = 50 MLs. The F overlayer with magnetization per- 
pendicular to the F|TI interface is assumed to induce surface 
gap Asurf in the helical state of TI via the magnetic proximity 
effect. 

ates between massless and massive Dirac fermions on the 
surface of the TI, it cannot be used to reliably detect the 
latter since it would diminish if bulk charge carriers are 
present within the TI slab. 

In the EM radiation detection schemes, the TI slab 
which is semi-infinite is described by topological field 
theory (TFT) [2], or alternatively one considers TI slab 
of finite thickness with its bulk and surface conductivi- 
ties pjj. The latter includes both the longitudinal and 
transverse (Hall) contributions [H^ where the role of 
{ae/47r^)E • B term (E and B are electromagnetic fields 
and a is the fine structure constant) in TFT is completely 
replaced by the appearance of explicit surface Hall con- 
ductivities (quantized or non-quantized) that influence 
the boundary conditions for EM wave propagating per- 
pendicularly through the slab [11]. Our device in Fig. [l] 
for measuring currents and voltages is analogous to this 
viewpoint where the presence of TI surface, with strong 



SOC and possibly opening of Ag^rf gap, in contact with 
the precessing F layer rectifies pumped spin current into 
charge current. The major advantage of our device is 
that due to perpendicular-to-plane transport geometry 
it can operate even if there are bulk charge carriers in 
the interior of the TI slab. 

We now explain details of our formalism to obtain 
pumped charge current and corresponding Vpump- The 
semi-MT J in Fig. [l] is modeled on a simple cubic lattice 
where monolayers (MLs) of different materials (N, F, TI) 
are assumed to be infinite in the transverse direct ion. 
The TI slab has thickness dxi = 5 MLs and the F layer 
has thickness dp = 50 MLs. The F and N layer are 
described by a tight-binding Hamiltonian with a single 
s-orbital per site 

n,crcr', k|| 

-7 XI (4a,k||Cn+l,<T,k|| +H.C.). (1) 
n,cr,k|| 

The operators cj^^ (cna) create (annihilate) electron with 
spin a on ML n with transverse momentum ky within the 
ML. The in-plane kinetic energy is £n,k|| = — 27(cos kya-\- 
cos/c^a), whose effect is equivalent to an increase in the 
on-site energy, and the inter-layer hopping is 7 = 1.0 eV. 
The coupling of itinerant electrons to collective mag- 
netization dynamics is described through the material- 
dependent exchange potential A^ = 0.2 eV (A^ = 
within semi-infinite ideal N leads), where & = {cFx^^yi ^z) 
is the vector of the Pauli matrices and [(jQ,]crcr' denotes 
the Pauli matrix elements. The Hamiltonian in Eq. ([T]) 
is time-dependent since the unit vector m(t) along the 
local magnetization direction within the F layer, which 
is assumed to be spatially uniform, is steadily precess- 
ing around the z-axis with a constant cone angle 9 as 
illustrated in Fig. [T] 

The minimal model for the slab of 3D TI, such as 
Bi2Se3, is the effective tight-binding Hamiltonian with 
four orbitals per site [17]: 



Hi' 



E 

n,k|| 



-^n,k|| 



B. 



2a 



^3 Cn+l,k|| + H.c. 



n,k|| 



CI + (i(k|| )ro H (Fi sin kxa + T2 sin kya) ^^m(t) • & 

(2) 



^n,k|| J* 



It yields the correct gap size in the bulk and surface dis- 
persion while reducing to the continuum k • p Hamilto- 
nian in the small k limit. Here c = (c+^, c+|, c_^, c_|)-^ 
annihilates electron in different orbitals, <i(k)|| = M — 
2B/a^ + 25(cos/c^a + cos/c^a-2)/a^ Ti {i = 0,1,2,3) 



are 4x4 Dirac matrices and 1 is the unit matrix of the 
same size. The numerical values of parameters are cho- 
sen as: M = 0.3 eV; A = 0.5 aeV; and B = 0.25 a^eV. 
The Fermi energy of the whole device is set at either 
Ef = 3.01 eV or Ef = 3.35 eV in order to force it to lie 
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inside (which is topological transport regime) or outside 
of A^^f , and the bottom of the band of the TI slab is 
shifted by C = 3.0 eV. 

Figure [sja) shows the local density of states (LDOS) 
on the TI surface within the semi-MTJ where the energy- 
momentum relation bears little resemble to the Dirac 
cone for strong coupling between the F layer and TI 
determined by the corresponding inter-layer hopping 
7c = 0.4 eV (or higher). This is due to flooding of the 
F|TI interface by evanescent wavefunctions which origi- 
nate from the F layer and penetrate into the bulk gap of 
TI while exponentially decaying in space, which is a coun- 
terpart of the well-known metal induced gap states in 
metal-semiconductor Schottky junctions. In other words, 
for strong F to TI coupling the TI surface state hybridizes 
with F metallic states so that sharpness of the Dirac cone 
is lost [H]. 

To avoid such situation, we assume smaller hopping 
7c = 0.15 eV where the corresponding LDOS is shown in 
Fig. |3|b). The weak F to TI coupling can be achieved 
by growing an ultrathin layer of a conventional band in- 
sulator, such as In2Se3 with large bandgap and good 
chemical and structural compatibility with Bi2Se3 where 
sharp heterointerfaces have already been demonstrated 
by molecular-beam epitaxy growth [19]. This dielectric 
layer will also affects the value of the proximity induced 
surface energy gap Agurf- In addition. Fig. |3] shows that 
gapped surface state of TI can penetrate into the bulk 
of TI as evanescent wavefunction decaying over the first 
three MLs which effectively dopes the bulk. 

Although the theory of spin pumping by F|N [13j or 
F|I [ini [20] interfaces has been developed to a great ex- 
tent, it is only very recently that strong SOC effects 
at the interfaces of F layers have been taken into ac- 
count using scattering approaches |21J or time-dependent 
NEGF [14, 22^. Since pumped charge and all components 
of pumped spin current tensor are time-dependent in the 
presence of SOC NEGF formalism is more advan- 
tageous choice because it gives from the outset experi- 
mentally measurable current averaged over one period. 
In contrast, scattering approaches [21 require to com- 
pute current at all times (or, in practice, over a discrete 
time grid) during one period of microwave oscillations 
and then find its average, which can be difficult in tun- 
neling structures where current amplitude can be several 
orders of magnitude larger than its average value. 

To understand charge pumping in the junction de- 
picted in Fig. [l] we employ our recently derived fl^ ex- 
act multiphoton solution to double-Fourier-transformed 
NEGFs in the presence of external time periodic fields. 
Using the operators cj^^ {cna) which create (annihilate) 
electron with spin a on ML n, we can introduce the two 
fundamental objects [23 of the NEGF formalism — the re- 
tarded G^^ {t,t') = -ie{t - n{{cna{t),cl^^{t')}) and 
the lesser G^^/ {t,t') = i{cl^.{t')cna{t)) GF that de- 



scribe the density of available quantum states and how 
electrons occupy those states, respectively. Here (. . .) 
denotes the nonequihbrium statistical average [23] . 

When the device Hamiltonian depends on time ex- 
plicitly, as is the case of + H^^ above, one has to 
work with both times. To solve the equation of mo- 
tion for G^'^r {t^t') cind the Keldysh integral equation 
for G^^T^ (^7^0 is advantageous to switch to a more 
convenient representation {h = 1) |24j 

+ 00 +00 

G''<{t,t')= J ^ j |Je-'^*+'^'*'G'-'<(^,E')- (3) 

— OO — oo 

Due to the Floquet theorem, the double-time Fourier 
transformed retarded GF, G'^{E^ E')^ must take the form 
G^(£;, E') = G^(£;, E + Nlj) = G^(£;). The coupling 
of energies E and E + Nuj^ where TV is integer, signi- 
fies how multiphoton processes assist the generation of 
pumped current. Since the junction in Fig. [l]is transla- 
tionally invariant in the transverse direction, we addition- 
ally perform spatial Fourier transform in the yz-plane, so 
that multiphoton retarded GF depends on ky . 

Assuming absence of electron interactions with 
phonons, magnons, or other electrons, solving [El — 
Hk|| - tl^^ {E + r2)]Gk|| {E) = i for the multiphoton re- 
tarded GF is sufficient to express the time-averaged total 
pumped charge current in N lead p = L, as: 

Ip = ^/rfk||Tr{f^,k„[17,G;„(i?^)fk„(i?f)]}, 

[n, Gk„ (i?F)f k„ (Ef)] = f2Gk„ (i?F)f k„ (Ef) 
-Gk„(^F)fk„(^F)f2, (4) 

where adiabatic limit u <C Ep satisfied by microwaves 
is taken into account. The symbol A means that ma- 
trix acts in the Hilbert space Hei Hph, where the 
dimension of the Hilbert space of photons Hph is infi- 
nite [14] and the dimension of Hei is equal to the num- 
ber of F and TI MLs, dp + <^ti- The unit matrix in 
the Hilbert space of a single electron Hei is 1, so that 
the unit matrix in Hei (8) Hph is denoted by 1. Here 
Ct = diag (• • • , — 2cj1, — cjI, 0, cjI, — 2co'1, • • • ) is the diag- 
onal matrix, T = Tp is the level broadening matrix 
due to the coupling of the F|TI bilayer to N leads. Since 
the trace in the integrand, Tr = TreiTrph, is summing 
over contributions from different photon exchange pro- 
cesses, the denominator includes 2A/'ph to avoid double 
counting. Note that the part of the trace operating in 
Hph space ensures the current conservation in our so- 
lution to NEGF equations. We find that A^ph = 10 is 
sufficient to obtain asymptotic value of pumped charge 
current. The pumped charge current is then trans- 
lated into the dc pumping voltage in an open circuit 
via Vpump = Ip/G{9) which is the quantity measured 
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in recent experiments on F|I|N semi-MTJs driven by mi- 
crowaves [15 . 

Note added. During the preparation of this manuscript 
we noticed Ref. ^5] where charge pumping was found for 
the same device as in Fig. [l] but with current flowing on 
the surface of TI (i.e., within the ?/2;-plane) rather than 
perpendicularly through such interface. Another major 
difference is that TFT treatment in Ref. [25 , requires per- 
fectly insulating bulk, while our vertical junction setup 
can operate even if there are excessive bulk charge carries 
in the interior of the TI slab. 
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